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Abstract 

We consider a small and fixed number of fermions (bosons) in a trap. The 
ground state of the system is defined at T = 0. For a given excitation energy, 
there are several ways of exciting the particles from this ground state. We 
formulate a method for calculating the number fluctuation in the ground state 
using microcanonical counting, and implement it for small systems of nonin- 
teracting fermions as well as bosons in harmonic confinement. This exact 
calculation for fluctuation, when compared with canonical ensemble averag- 
ing, gives considerably different results, specially for fermions. This difference 
is expected to persist at low excitation even when the fermion number in the 
trap is large. 
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I. INTRODUCTION 



The traditional approach to determining the number fluctuations in a given quantum 
state relies on the statistical description of the system based on the grand canonical ensem- 
ble (GCE). For bosons, if < Ni > is the average occupancy of a single orbital i, then the 
mean-square fluctuation in GCE is given by JIl < (AA^j)^ >=< Ni > (1+ < Ni >). As the 
temperature T — > for the boson system, however, there is macroscopic occupancy in the 
ground state, with the result that the relative fluctuation < (AA^q)^ >^^^ /N 1. This is 
clearly absurd, since at T = 0, none of the bosons (in the noninteracting case) should be 
in the excited orbitals, so the fluctuation should vanish. This problem, flrst discussed by 
Fujiwara et al. 0, is of topical interest with the discovery of Bose-Einstein condensation 
in a dilute gas of alkali atoms conflned in a trap The problem of number fluctuation 
from the ground state of a system of conflned fermions is also of interest due to the recent 
developments in the fleld of quantum dots [Q. The fluctuation in Bose-Einstein systems 
has been addressed in a series of papers both in micro-canonical and canonical ensem- 
ble framework. A clear and unambiguous deflnition of the fluctuation has been given for 
calculating the ground state occupation number fluctuation by Grossman and Holthaus 
through the "microcanonical entropy" of an N-particle Bose system. Detailed comparisons 
between the two methods for noninteracting bosons in a harmonic trap have been made by 
Holthaus and Kalinowski 0, where analytical expressions for the two are also given using a 
modifled saddle-point approximation. 

Our focus in this paper is on the exact calculation of the ground state fluctuation for a 
small number of fermions conflned in a trap. For comparison, the corresponding calculations 
are also done for bosons. To this end, a combinatorial method based on microcanonical 
counting is developed in sect. (II), and calculations are made for particles in one and two- 
dimensional harmonic traps. These exact results are then compared with the corresponding 
canonical ensemble calculations. The two methods yield substantially different, though 
qualitatively similar results. As explained later, we expect this deflciency of the canonical 
ensemble averaging method to persist for fermions at low excitation, even when the fermion 
number in the trap is large. In sect. (Ill) we show that even though the canonical entropies 
for noninteracting bosons and fermions in a one-dimensional harmonic trap are identical, 
the number fluctuations in the ground state are vastly different. The numerical results are 
discussed in sect. (IV). 

II. FLUCTUATIONS IN THE INDEPENDENT PARTICLE MODEL 

We assume an independent particle model where the one particle partition function is 
given analytically, or may be computed. We flrst give the well-known recipe of calculating 
the A^- particle quantum canonical partition system from this. At zero temperature all the 
particles are in the ground state. At a flnite temperature, however, the excitation energy may 
be shared in many different ways amongst the particles, with the result that the population 
in the original ground state is not flxed, although the total number, A^, is still the same. 
Our objective here is to deflne and calculate this fluctuation in the ground state occupation 
as a function of excitation energy or, equivalently, temperature. 
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The canonical partition function for bosons and fermions in any space dimension may 
be written as 



zr = {±f E n^^^^^^rr^' (1) 



N 



71 ■ ' 

ni,n2,.--,'T.jv i=l J' 



where Zi{(3) is the single particle partition function, (3 is the inverse temperature and the 
upper and lower signs refer to bosons and fermions respectively. The sum over the set of 
integers is constrained by the relation 

N 

E =N. (2) 

i=i 

The above formula allows us to write a general recursion relationship for the canonical 
partition function 

1 ^ 

Zn'^ = -^T.i±r-''Z,{nP)Z^f^{P) (3) 

n=l 

for bosons (+) and fermions (— ). We note that in the above recursion relation Zo is formally 
taken to be unity for consistency. 

In order to perform explicit calculations, we specialise to the case of a harmonic oscillator 
in d-dimensions. The single particle partition function is given by, 

= 1^1^ (4) 

where 

X = exp{— (3Tiuj) . (5) 

The canonical partition function for a system with particles is then computed using 
Eq.(|I]) and is given by 

ZN = x^'^l'P^{x)f{^^^^^^ (6) 

where Pn{x) is a polynomial in x which depends on the dimension and the statistics of the 
system. We shall use the notation Pn{x) = Bj^{x) and Pn{x) = F]y{x) for bosons and 
Fermions where necessary. The polynomial may be calculated using the recursion relation 
in Eq.(|) 

The recursion relation above should be used with the condition Po{x) = 1 for both bosons 
and fermions. We further note that in one-dimension Bn{x) = 1 for bosons and Fn{x) = 
^N{N-i)/2 fermions. They are however more complicated in higher dimensions ||T3| . 
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A. Fluctuations from Microcanonical Counting 



We first define the fluctuation in particle number from the ground state at a given 
excitation energy through a set of counting rules. Again we first write down the general 
formulae for a given a set of discrete energy levels, and then specialize to the harmonic trap. 
The single particle partition function may be written as 



oo 



Zi(/5) = x^o^x^% (8) 
i=i 

where x = and ej,j = 0, .., oo are the single particle energies. It is understood that (3 in 
the exponent defining x has been multiplied by a characteristic energy scale of the system, 
and similarly ej has been divided by the same, which we put to unity for convenience. For 
the harmonic oscillator, this energy scale is hu, and x is given by Eq.(||). Substituting this 
into Eq.(0), and expressing Zn in a power series in x, we obtain 

oo 

Zj, = x^'^j:n{Et\N)x''t\ (9) 
fe=i 

where the A^-particle eigen-energies = Eq + E^^^ form an ordered set, with Eq and E^^^ 
denoting the ground state energy and the excitation energy with respect to the ground state 
respectively. The expansion coefficient VL{E^j^^\n) denotes the number of possible ways of 
distributing the excitation energy E^^"* in utmost N particles 
Furthermore we may write fi(4''^\Ar) as 

N 

mt\N)= Y: io{Et\N^.,N) (10) 

where uj{Ej^^\ Nex, N) denotes the number of possible ways of distributing the excitation 
energy Ej^^"^ amongst exactly Nex particles. Hence the probability of exciting exactly Nex 
particles from an A^-particle system at an excitation energy E^^'' is given by 

.(e.) ,n _ ^ {Et^\Nex,N ) 

n{Et\N) 



piEt\N^x,N) = '^'2Jie7ZJ Nex = 0,...,N. (11) 



By definition this probability is properly normalised. Further the probability has the fol- 
lowing properties: 

p{0,N,x,N)=6oN.^, (12) 
piEt\N,x,N) = N,x>N. (13) 

The number fluctuation in the ground state of the system may now be defined in terms 
of the moments of the probability distribution given above. We first define the moments: 

N 

<N,x> = E N,xPiEt\N,x,N), (14) 

N 

<Nlx> = E Nlp{Et\N^x.Nl (15) 
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and the number fluctuation from the ground state is given by, 



= <7V2>_<7Vo>2 , (16) 

since < iVo > + < N^x >= N. 

Few remarks are in order here: The above definitions apply equally well to bosonic and 
fermionic systems. The fluctuation in the number of particles from the ground state is ex- 
pressed here as a function of the excitation energy with respect to the ground state. In the 
case of bosons this is just the fluctuation from the lowest energy single-particle state, where 
as for fermions it is the number fluctuation across the (zero temperature) Fermi energy. For- 
mally the above expressions complete the necessary basic definitions for further analysis. In 
actual practice the counting of the number of possibilities to which a system may be excited 
is non-trivial. In the following, we give an example of bosons and fermions in a harmonic 
oscillator potential for the exphcit calculation of uj(eI^^\ Nex, N). Henceforth, we consider 
the excitation energy from the ground state to be nhcu, and denote the corresponding mi- 
crocanonical multiplicity to be a;(n, N^^, N). 

(a) Bosons in a d-dimensional harmonic oscillator. 
In the case of a harmonic oscillator trap, all the single particle energy levels in a given 
shell have the same energy. The shell may be characterised by the index k with energy 
{k — 1 + d/2)tuj, k = 1,2, ..oo, and degeneracy Qk- When the excitation energy is zero 
(equivalently at zero temperature) all the particles are in the ground state. Now consider 
exciting Nex particles from this ground state that share n quanta of energy. Let there be mi 
bosons with 1 quantum of excitation energy, m2 bosons with 2 quanta of excitation, and so 
on. We then have 

n^mi + 2m2 + ... + {n+1 - Nex)mN^^ 

and 

Nex = mi + m2 + ... + mN,^. 

This is just a way of partitioning an integer n into exactly N^x partitions with denoting 
the number of times an integer i appears in a partition ( it is taken to be zero if some integer 
between 1 and n + 1 — N^x does not appear in the partition). All the bosons occupy 
the state i since i — 1 is the ground state. The number of ways these rrii bosons are 
distributed in the state i -|- 1 is given by the counting rule 

(g,+l+mi-l)^^^_ (17) 

The set {rrii) denotes the set of all partitions of n for a given N^x- The micro canonical 
distribution uj is then given by 

Nex 

uj{n, Nex, iV) = E n ^''^'^""'-'^Cm,. (18) 

{mfc} k=l 

The formula for uj is especially simple in one dimension since there is no degeneracy, Qk — ^ 
for all k. We have. 
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u{n,N,,,N) = E 1, 



(19) 



and the number of possibilities are determined by the number of ways an integer n can be 
partitioned into N^x integers, that is the number of ways we can write, 

n = m + ^2 + ria + ... + tln,^, ni > n2 > ...riN^^. 

Once the tu's are known the probabihty distribution may be calculated using the Eq.(|ll]) 
and hence the fluctuations as a function of the excitation energy. In order to express this 
as a function of temperature, one may use the relation between excitation energy above 
the ground state and the temperature to be discussed presently (see Eg. (p3D . We, however, 
prefer to present all results as a function of the excitation energy only. Since N^^ < N 
always, we can give an alternative expression for as 

uj{n, Nex, N) = n{n, N^x) - ^{n, N^x - 1). (20) 

While this allows an alternative and probably a simpler way of calculating probability dis- 
tribution, it is valid only for bosons (see the discussion following Eq.(p4[) given later). 

(b) Fermions in a harmonic-oscillator (closed shell). 
For simplicity, the fermions are taken to be spinless, and filling an integral number of shells 
in the ground state. The determination of microcanonical distribution uj{n, N^x, N) is more 
involved than in the case of bosons. The combinatoric rules depend on the space dimension. 
Using the same notation for the harmonic oscillator single particle energies in a given shell, 
we calculate the number of possibilities by which Nex particles can be excited from the 
ground state. Note that removing fermions from the ground state leaves as many holes 
in ground state. As a result the u depends not only on the distribution of fermions in 
the excited states (as in the case of bosons), but also on how the holes are distributed in 
the ground state. For simplicity we choose a completely closed shell system, the rules may 
be appropriately modified for open shells. We assume that N-fermions form a closed-shell 
system at T = with shells up to k (corresponding to fermi energy Ep = (A; — 1 -|- d/2)%uj) 
filled. Thus Nf,x holes are distributed in states from i = 1 xvp to i = k. Let hi denote the 
number of holes in the state i such that 

Nex = hi + h2 + ... + hk. 

The number of ways hi holes may be created in the zth state is given by ^^Ch-. Hence the 
number of ways in which Nex holes may be created in the ground state is given by the 
product over the number of ways in which the holes may be distributed in the k orbitals, 
i.e.. 

Now consider exciting Nex particles from this ground state sharing n quanta of energy. An 
allowed configuration is one in which each and every one of these Nex particles is found in 
states above the fermi energy, with the shell indices ranging from [k + 1) up to (A; + ra), such 
that their excitation energies add up to yield the total Eex- This complicates the counting 
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rules for fermions as compared to bosons. We shall denote the occupancy of orbitals for the 
excited particles by rrij, where i = k + 1, k + n. The number of ways the femions are 
distributed in the state k + i is then given by the counting rule 

i9k+i)r< 

The micro canonical distribution u is then given by 

k k+n 

a;(n,iVe,,iV) = E E 11 ^''^^.^ U ^''^'^Cm., (21) 

{mi}{hj}j=l i=k+l 

where Ngx = J2i and the microcanonical multiplicity uj is obtained by summing over all 
the allowed possibilities such that the sum total of the excited quanta is exactly n. Once the 
w's are known the probability distribution may be calculated using the Eq.(^) and hence 
the fluctuation as a function of the excitation energy. 

We compare these results with the fluctuations obtained by the canonical ensemble av- 



eraging method of Parvan et al. as detailed below. It is our objective to see how close 



are the results for the ground state fluctuations calculated by the two methods. 



B. Fluctuations from canonical ensemble averaging 

In statistical thermodynamics, a macro-state at a given energy may be built up by many 
combinatorial ways from the micro-states, and this is denoted by the term multiplicity of the 
macro-state. As we saw in the preceding sect. (II A), the multiplicity VL{EI^, N) was defined 
by Eq.(|lQ|) through this counting method. In the canonical ensemble, we may alternately 
define 

fi(T, A^) = exp[5(T, A^)] = exp[(t/ - F)/T] = Z^{x) . (22) 

In the above, the internal energy t/(T) is determined as usual from the canonical partition 
function, and excitation energy at any temperature is given by 

^(ex) _ ^(^^ _ ^(Q) _ (23) 

We may therefore compare the calculated quantities from the canonical and the microcanon- 
ical ensembles as a function of the excitation energy. 

Comparing J1(T, A^) given by Eq.(p^ from the canonical ensemble with the series (|]), we 
see that it is as if only one term from this series is picked in the ensemble averaging. This 
is realized for a large number of particles, since the multiplicity fi(i?f^, A^) increases rapidly 
with the excitation energy, where as the factor x^^ decreases exponentially. In this paper, 
we focus on systems where A^ is not large, specially for fermions. It is therefore interesting 
to examine the differences in the results of the calculation made by the two methods. 

Once the multiplicity is obtained, the fluctuations may be defined through a temper- 
ature dependent probability distribution 

Mr.iv.,.iV)^ "'^-^"'^^"^^-^"-') . iV.^O,...,iV. (24) 
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Note that ft{T, N^x) is determined by the canonical partition function of Nex particles by 
replacing by N^x in Eq.(^). For bosons, the ground state of every such system consists 
of a single level, and the definition given by Eq.(p^ is correct. For fermions, however, the 
Fermi energy of a system defined by Nf.x particles is less than the Fermi energy of the full 
system with particles. The fluctuations are defined with respect to the ground state with 
the Fermi energy corresponding to the full system. Therefore Eq.(p4D cannot be applied to 
fermions. A consistent definition of ground state fluctuations applicable to both fermions 
and bosons, and which coincides for the bosonic calculation with the probability definition 
(P^), may however be given by the ensemble averaging method |rT,T2|. We summarise the 
method below. 

The canonical partition function in Eq.(P) may be written in the occupation number 



representation as 12 



^7v=Ena^''"N (25) 

{nfe} k 

where we have used the fact that the energy of the A^- particle system for a given set of 
occupancies {n^} is given hy E = Y^k^k^k- The occupancy = 0, 1 for fermions and may 
take any value up to A^ for bosons. At finite temperatures, using the recursion relation in 
Eq.(^, and some nontrivial algebra, the ensemble-averaged moments of the occupancy 
may be written as [|12| 



1 ^ 



<nk> = ^ Y^i^y^'^'^'Z^.^, (26) 
Zn j=i 

1 ^ 
Zn 

N N-j 

j=l 1=1 



Zn 



where the upper and lower signs refer to bosons and fermions respectively. 
The ground state fluctuation is now given by, 

5Nl = Y.{<nl>-<n,>'). (28) 

k 

The sum runs through all the allowed k values in the ground state defined at zero tempera- 
ture. Note that the moments of occupancy are defined without any reference to the ground 
state. The fluctuation is then obtained by summing over the single particle states occupied 
at zero temperature. In the counting method as outlined earlier, the fluctuations necessarily 
refer to the ground state. 

Applying the above equations to a system of bosons in a harmonic trap, we get 

1 ^ 

<no> = —Y.x'''/^Zn-,, (29) 
Zn ,=i 

1 ^ 

< ^0 > = ^ E(2J - l)x^''''ZN-y (30) 
Zn 
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The ground state fluctuation is given by Eg. (pH]) , 

6N^ = {<nl>-<no>^). (31) 

Similarly, the thermodynamic expression for the ground state fluctuations for the 
fermionic system in a harmonic oscillator is obtained from Eq.(p6D and Eg. (p?]) : 

<ns> =^ j:{-iy^'9sX^^'-'^'^'^Z^.„ (32) 
Zn 

1 ^ 

<nl> =^ Ei-^y^'iJds - U - l)g!)x^^'-'^'/'^ZN-„ (33) 
Zn 

where Qs refers to the degeneracy factor of the the state s. For the N-fermion system 
that forms a closed shell system at T = with fllled shells up to k, with Fermi energy 
Ep = (A; — 1 + d/2)hu!, the ground state fluctuation is given by Eg . (|28|) : 

k 

^iVo =E {<<>-<n,>'). (34) 

Before using these formulae for the calculation of ground state occupancy and fluc- 
tuations, we discuss below the specially interesting case of the one-dimensional harmonic 
confinement where the canonical and the microcanonical ensembles yield vastly different 
results. 



III. FLUCTUATIONS IN A ONE-DIMENSIONAL HARMONIC TRAP 

The one-dimensional Harmonic trap is specially interesting because even though the 
canonical entropies for bosons and fermions are identical, the number fiuctuations from 
the ground state are very different. This may be seen by writing the canonical iV-particle 
partition function as 

Z^ = x^^(^-'y'^^/'f[-±-, (35) 
j=i U ~ ^ J 

with g = for bosons, and g = 1 for fermions. Actually, for other positive values of g, the 



above form is the exact partition function for the so-called Calogero- Sutherland model [10 



where the N particles interact pair-wise by a potential ^J2i<j{xi — For bosons, the 

dimensionless parameter g is in the range < g < 1/2, while for fermions g > 1/2. The 
special values g = 1{0) give noninteracting fermions (bosons). The effect of interaction has 
only been to shift the energy of every state by the same amount, which is absorbed in the 
prefactor. It follows from Egs.(|^,^) that we may write 

oo 

Zj, = :,9N(N-l)/2+N/2 g ^^^^ ^ (gg) 

n=0 

and that f2(r;,, A^) is independent of the parameter and is the same for bosons and fermions. 
Since it is the logarithm of VL{n, N) that determines the canonical entropy of the system at 
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an excitation energy of n quanta, it follows that the entropy is independent of g. The same 
result is true if one calculates the ensemble averaged entropy at a given temperature. This 
may be easily verified by using the relation F = —In Zj^j / (3 for the free energy, and then 
calculating the entropy S' = — 

For the microcanonical calculation of the fluctuation, we need to calculate the micro- 
canonical multiplicity uj{n, Nex, N), which is the number of ways of distributing the n exci- 
tation quanta amongst exactly Nex particles. Although the relation 



is obeyed both by fermions and bosons, and the LHS of the above equation is the same 
for both, the microcanonical counting of cu's are very different for the two cases. In the 
appendix, this is illustrated explicitly for = 3. Thus the fluctuations for the bosonic 
and fermionic cases differ substantially when the exact counting method is used. As is to 
be expected, the fermionic fluctuation is considerably suppressed compared to the bosonic 
system. The same qualitative conclusion may also be reached by using Eqs.(|3y, ^Sj) based 
on canonical ensemble averaging rather than exact counting. In the next section, we display 
these results, as well as the results for the two-dimensional harmonic oscillator. 



In Fig.(l), we display (a) the ground state occupancy < Aq > and (b) the relative fluc- 
tuation < 6Nq > /N for A^ = 15 noninteracting particles in a one-dimensional harmonic 
oscillator potential. The canonical and the exact microcanonical results are compared as a 
function of the excitation energy. The canonical method gives results in close agreement with 
counting for the ground state occupancy < Nq >, but overestimates the relative fluctuation 
substantially. In one dimension, the number of microcanonical possibilities uj{n, N^x, N) is 
very restricted at low excitations due to the non- degeneracy in the single-particle energy 
levels and the Pauli exclusion principle. This results in (i)< A'q > for fermions getting 
depleted more slowly than bosons, and (ii) reduced fluctuation. In Fig. (2), the same quanti- 
ties are displayed for particles in a two-dimensional harmonic oscillator potential. It is well 
known from previous work that there is a peak in the relative fluctuation somewhat 
below the critical temperature for bosons. For fermions, such peaking is absent, specially 
in the exact calculation. We thus see that the presence or absence of a peak in the number 
fluctuation of the ground state may signal phase transition, or its absence. 

The microcanonical method of exact counting for fermions is computationally very time- 
consuming because, unlike bosons, Eq . (^4|) cannot be used. We have therefore restricted 
the fermionic calculations to only up to A^ = 15 particles. As we see from Figs (1-2), 
there is considerable difference in the results for the relative fluctuations for small particle 
numbers. For fermions, we expect this difference to persist at low excitations even when A^ 
is large. This is because at low excitations, only a small fraction of the fermions near the 
Fermi sea can be excited, so the effective number of fermions contributing to the number 
fluctuation remains small even when the system is large. For bosons, there is no difficulty 
in performing the microcanonical calculation for a large number of particles. In Fig. (3), 



N 




(37) 



IV. RESULTS AND DISCUSSION 
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we show the bosonic occupancy and the relative fluctuation as a function of the number of 
excitation quanta for = 100. We note that the occupancy < Nq > from the canonical and 
exact microcanonical methods is practically identical, although the ground state fluctuation 
continue to differ. 

We have made the comparison between the microcanonical and canonical results as 
functions of the excitation energy, rather than temperature. In the microcanonical method, 
the counting of the possibilities is done for a given excitation energy, the latter being the 
natural variable. In the canonical formalism, of course, a mapping from the excitation 
energy to temperature can be done using Eg. . We may also note that starting from the 
canonical Eq.(p2D and using Eq.(p4D for bosons, one can obtain the analytical expression for 
the low temperature dependence of the fluctuation, as, for example, given by Eq.(14) of p. 
It remains a challenging problem to obtain similar microcanonical relations for fermions. 

This research was supported by the Natural Sciences and Engineering Research Council 
(NSERC) of Canada. M.T. would like to acknowledge receiving an NSERC post-graduate 
scholarship, and M.V.N.M. enjoyed the hospitality of the department of Physics and As- 
tronomy of Mcmaster University, where the work was done. 
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FIGURES 



FIG. 1. (a) Plots of the ground state occupancy < Nq > versus the excitation energy E^x for 
N = 15 bosons (fermions) in a one-dimensional harmonic oscillator. E^.c is given in units of the 
oscillator spacing. The microcanonical results are displayed according to the legends in the inset 
box. (b) Plots of the relative ground state fluctuation < 6Nq > /N as a function of E^x for the 
same systems as in (a). 

FIG. 2. (a) Plots of the ground state occupancy < Nq > versus the excitation energy E^x for 
A'^ = 10 bosons (fermions) in a two-dimensional harmonic oscillator, according to the legends in 
the inset box. (b) Plots of the relative ground state fluctuation < SNq > /N as a function of E^x 
for the same systems as in (a). 

FIG. 3. (a) Plots of the ground state occupancy < Nq > versus the excitation energy Eex for 
N = 100 bosons in a one-dimensional harmonic oscillator, (b) Plots of the relative ground state 
fluctuation < SNo > /N as a function of E^x for the same system as in (a). 



12 



V. APPENDIX 



Calculation of uj{n, N^.^, N) in one-dimensional harmonic oscillator, 
(a) Bosons 

For illustration, in Table 1, we display the simple case of = 3 bosons, with the the number 
of excitation quanta n < 7. Since there is no degeneracy, the combinatorial factor given by 
(p!7|) is unity. Therefore ci;(n, N^^, N) in this case is just the number of distinct ways in which 
the integer n may be partitioned amongst exactly Ngr^ identical bosons (A^^e^: < 3). In Table 
1, n increases from left to right, and increasing values of N^^ are tabulated in a vertical col- 
umn. The integer in each box is the corresponding u!{n, N^^, N), with the distinct partitions 
of n listed in brackets below it. For example, in the box under {n = 4, N^x = 2), we see that 
a;(4, 2, 3) = 2, and the two distinct partitions of 4 are (3 + 1) and (2 -|- 2). In the last row is 
listed Q{n, N), that is obtained by adding the uj{n, N^x, A^)'s in each vertical column. Note 
that these check with the coefficients in the expansion of the 3-boson canonical partition 
function : 

Z^ = l + x + 2x^ + 3x^ + 4x^ + + 7x^ + + .... 

(b) Fermions 

Same as for bosons, except that two more vertical columns have been added under the 
headings Li and n^m- For three spinless fermions, the lowest three energy levels (1, 2, 3) (in 
increasing order of energy) are occupied at T = 0. The column under Li lists the possible 
configuration of holes in these levels for an excitation energy of n > n^m quanta, the latter 
being listed in the adjacent column. For example, for Nex = 2, (3i2i) under the column 
Li denotes a two-hole configuration, with one hole in level 3, and another in level 2. The 
minimum energy required for this is n^in = 4 in units of hu. 
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TABLES 



TABLE L Tabulation of bosonic io{n, N^x, N) for N=3 



n= 


1 


2 


3 


4 


5 


6 


7 




1 


1 


1 


1 


1 


1 


1 


Lo{n,Nex,N): Nex=2 





1 

(1+1) 


1 

(2+1) 


2 

(3+1) 
(2+2) 


2 

(4+1) 
(3+2) 


3 

(5+1) 
(4+2) 
(3+3) 


3 

(6+1) 
(5+2) 
(4+3) 


Nex=3 








1 

(1+1+1) 


1 

(2+1+1) 


2 

(3+1+1) 
(2+2+1) 


3 

(4+1+1) 
(3+2+1) 
(2+2+2) 


4 

(5+1+1) 
(4+2+1) 
(3+3+1) 
(3+2+2) 


n{n, N) = 


1 


2 


3 


4 


5 


7 


8 



TABLE XL Tabulation of fermionic w(n, Nex, N) for N=3 





n 


Li 




1 


2 


3 


4 


5 


6 


7 


3i 


1 


1 


1 


1 


1 


1 


1 


1 


Nex = l- 2i 
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